Fermi-liquid, non-Fermi-liquid, and Mott phases in iron pnictides and cuprates 
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The role of Coulomb correlations in the iron pnictide LaFeAsO is studied by generalizing exact 
diagonalization dynamical mean field theory to five orbitals. For rotationally invariant Hund's 
rule coupling a transition from a paramagnetic Fermi-liquid phase to a non-Fermi-liquid metallic 
phase exhibiting frozen moments is found at moderate Coulomb energies. For Ising-like exchange, 
this transition occurs at a considerably lower critical Coulomb energy. The correlation-induced 
scattering rate as a function of doping relative to half-filling, i.e., 8 = n/5 — 1, where n = 6 for the 
undoped material, is shown to be qualitatively similar to the one in the two-dimensional single-band 
Hubbard model which is commonly used to study strong correlations in high-T c cuprates. In this 
scenario, the parent Mott insulator of LaFeAsO is the half-filled n = 5 limit, while the undoped 
n — 6 material corresponds to the critical doping region S c w 0.2 in the cuprates, on the verge 
between the Fermi-liquid phase of the overdoped region and the non-Fermi-liquid pseudogap phase 
in the underdoped region. 

PACS. 71. 20. Be Transition metals and alloys - 71.27+a Strongly correlated electron systems 
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I. INTRODUCTION 

The discovery of superconductivity in the iron-based 
pnictides 1 - — has recently led to an intense debate on the 
question whether Coulomb correlations in these mate- 
rials play a role of similar importance as in the high-T c 
cuprates£~— Moreover, in contrast to the effective single- 
band character of the 3d x 2_ y 2 states in the cuprates, the 
number of relevant d orbitals in the pnictides has also 
been a subject of discussion. In addition, because of the 
multi-band character of the Fe 3c? bands, the interplay of 
Coulomb repulsion and Hund's rule coupling should be of 
crucial importance for the pnictides. Finally, while it is 
generally accepted that the non-Fermi-liquid pseudogap 
properties in the underdoped regime of the cuprates are 
associated with the vicinity to the Mott insulating phase, 
in the pnictides it is not yet clear whether there exists a 
nearby Mott phase in the range of realistic Coulomb and 
exchange energies. 

To study the effect of Coulomb correlations in iron 
pnictides, various groups^^ ~ 18 i 20 ' 21 ' 25 ' 32 i 33 ' 37 have ap- 
plied single-site dynamical mean field theory 41 - (DMFT). 
Depending on the details of the single-particle Hamilto- 
nian and the magnitude of Coulomb and exchange ener- 
gies, weakly or strongly correlated solutions were found. 
In the present work we extend exact diagonalizatio n 42 i 43 
(ED) DMFT to five orbitals in order to investigate cor- 
relations effects in LaFeAsO as a function of Coulomb 
energy. The single-particle properties are described in 
terms of a five-band tight-binding Hamiltoniani 2 ^ Each 
d orbital hybridizes with two bath levels, giving 15 lev- 
els in total. The five baths are coupled indirectly via 
the interorbital 3c? Coulomb and exchange interactions. 
Thus, the Hilbert space is extremely large and finite-size 
effects are effectively reduced. The advantage of this 
multi-orbital ED/DMFT approach is that it is partic- 
ularly useful at low temperatures and that it can han- 
dle large Coulomb energies and full Hund exchange. As 



will be discussed below, the latter feature is of special 
importance in the pnictides since Hund's rule coupling 
leads to electronic properties that differ qualitatively 
from those obtained for the more approximate Ising-like 
exchange treatment. The sensitive role of exchange inter- 
actions in the pnictides was also noted in several previous 
papers 

The main result of this work is the identification of a 
paramagnetic Fermi-liquid to non- Fermi-liquid transition 
at moderate Coulomb energies, U c ~ 3 eV ( J — 0.75 eV), 
i.e., well below the overall width of the Fe 3c? bands, 
W ~ 4.5 eV. This incoherent metallic phase extends up 
to rather large values of U (> 6 eV if J is kept fixed 
at 0.75 eV). It is associated with the formation of lo- 
cal moments and with substantial low-frequency scat- 
tering rates in all 3c? bands. Below this transition, all 
bands exhibit strong correlation-induced effective mass 
enhancement. These properties are intimately related 
to the multi-band nature of LaFeAsO and the Hund's 
rule coupling among the Fe 3c? subbands. A similar spin 
freezing transition had been found recenly by Werner et 
alM: in a fully degenerate three-band model. Within the 
present five-band system, a Mott phase is not obtained 
below U — 6 eV. Thus, correlation effects in LaFeAsO 
appear to be related to the formation of local moments 
within the non-Fermi-liquid phase, and not to the vicin- 
ity of a Mott insulating phase. 

On the other hand, LaFeAsO readily turns into a 
Mott insulator at realistic Coulomb energies in the hy- 
pothetical limit of one-hole doping, i.e., n — 5. The 3g! 
bands then become half-filled and split into lower and up- 
per Hubbard bands. At intermediate hole doping, non- 
Fermi-liquid behavior dominates, while for electron dop- 
ing (n > 6) the system becomes a normal Fermi liquid. 
These results suggest a remarkable correspondence be- 
tween the multi-band compound LaFeAsO and the two- 
dimensional single-band Hubbard model. Indeed, if the 
correlation-induced scattering rate is plotted as a func- 
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tion of doping relative to half- filling, i.e., S = n/5 — 1 
(n = 6 for pure LaFeAsO), both systems exhibit the same 
sequence of phases for increasing <5: a Mott insulator at 
half-filling, a non-Fermi-liquid phase up to a critical dop- 
ing of the order of 5 C « 0.15 . . . 0.20, and a weakly corre- 
lated Fermi liquid for S > 5 C . In this scenario, the para- 
magnetic phase diagrams of iron pnictides and cuprates 
are strikingly similar, with LaFeAsO (6 = 0.2) located 
slightly on the overdoped side above critical doping. As 
a result, the system exhibits a clear asymmetry with re- 
spect to doping. Whereas electron doping beyond n = 6 
(5 > 0.2) reinforces Fermi-liquid properties, hole dop- 
ing n < 6 (6 < 0.2) enhances bad-metallic behavior. 
Although the real material is undoubtedly more com- 
plex because of doping-dependent single-particle proper- 
ties and the presence of antiferromagnetism^ at n = 6, 
we believe that the above picture nevertheless provides a 
useful new perspective for the role of correlation effects 
in iron pnictides in comparison with analogous physics in 
the cuprates. 

If exchange interactions among Fe 3d electrons are ap- 
proximated in terms of Ising-like exchange, i.e., by ne- 
glecting spin-flip and pair-exchange processes, the Fermi- 
liquid to non-Fermi-liquid transition still exists, but the 
critical value of U is shifted down to about 2 eV (assum- 
ing J = J7/4) and the transition is more abrupt. Similar 
qualitative changes from Hund to Ising exchange were 
found previously also for the Mott transition in two-band 
Hubbard models4^ 

Because of the strong hybridization between Fe 3c? and 
As 4p and O 2p states, there are indications that an 
accurate Wannier representation should encompass not 
only d but also p type basis functions, even if Coulomb 
correlations are explicitly only included among the 3d 
orbitalsi 18 ' 19 ' 21 ' 25 i 30 ' 32 A particular consequence of dp hy- 
bridization is that the effective 3d Coulomb interaction 
is considerably reduced. Moreover, Coulomb interac- 
tions among different d states are differently screened, 
giving rise to nonisotropic intraorbital and interorbital 
matrix elements . 19 ' 32 The present approach is general in 
the sense that these choices only affect the single-particle 
Hamiltonian and not the evaluation of the d electron self- 
energy matrix. Nevertheless, in this initial five-orbital 
ED/DMFT study we use, for simplicity, a purely d elec- 
tron tight-binding picture 2 ^ in order to elucidate the na- 
ture of the transition from Fermi-liquid to non-Fermi- 
liquid behavior. A more detailed investigation within a 
dp formulation is planned for future work. As in pre- 
vious papers,-^— £L£L2&2Zi2L21 we focus here on the 
paramagnetic phase. 

The outline of this paper is as follows. In Section II we 
discuss several theoretical details concerning the single- 
particle properties of LaFeAsO and the ED/DMFT pro- 
cedure that is used to evaluate the Fe 3d self-energy com- 
ponents. In Section III we analyze the results, with par- 
ticular focus on the differences obtained for Hund and 
Ising exchange. Subsection A discusses the neutral sys- 
tem, whereas the effect of doping and the analogy be- 



tween pnictides and cuprates are the subject of Subsec- 
tion B. Section IV contains the summary. 



II. MULTI-BAND ED/DMFT 

In this Section we briefly outline the theoretical de- 
tails of the multi-band ED/DMFT approach used in this 
work. The focus is on the role of Coulomb correla- 
tions within the Fe 3d subbands of LaFeAsO. The single- 
particle properties arc described in terms of the five-band 
tight-binding Hamiltonian H(k) which was recently de- 
rived by Graser et al^ for a single plane of Fe atoms 
from an accurate fit to the density functional results by 
Cao et al£L The low-energy part of these bands are in ex- 
cellent agreement with analogous calculations by Singh 
et a/.— The basis functions are d xz , d yz , d x 2_ y 2, d xy , 
and d 3z 2_ r 2, where the x, y axes point along Fe near- 
est neighbor directions. The first three of these orbitals 
comprise the t 2g subset, the remaining two orbitals rep- 
resent the e g subset. Hopping up to fifth nearest neigh- 
bors was included in the tight-binding fit. The onsite 
energies are E xz<yz = 0.13 eV, E x ?_ y i = —0.22 eV, 
E xy — 0.3 eV, and E z i — —0.211 eV. Thus the d xz ^ yz . xy 
levels lie 0.2... 0.4 eV above the d x 2_ y 2 z 2 levels. The 
hopping parameters are given in the Appendix of Rcf. 24 
and the one-electron band structure corresponding to this 
Hamiltonian is shown in Fig. 5 of this Reference. 

The Fe 3d density of states components are shown in 
Fig. 1. For symmetry reasons, the d xz and d yz compo- 
nents are degenerate. The widths of the t 2g and e g bands 
are approximately 3.0 eV and 4.0 eV, respectively. The 
total band width is about 4.5 eV. All Fe 3d bands ex- 
hibit a pronounced bonding-antibonding splitting, with 
a deep pseudogap at small positive energies, due to the 
hybridization with the neighboring LaAsO layers. In the 
absence of correlations, the occupancies of these bands 
are: n xz%yz = 0.58, n x 2_ y 2 = 0.53, n xy = 0.52, and 
n z 2 = 0.78. Note that these occupancies do not reflect 
the crystal field splitting among the onsite energies be- 
cause of the complex shape of the density of states com- 
ponents. In particular, the bands of d z 2 character are 
considerably more occupied than those of d x 2_ y 2 char- 
acter, despite the fact that E z 2 > E x i_ y i. A similar 
situation exists in the layer compound Naa;Co02, where 
n e < 9 > n ag although E e , g > E ag 42 

Previously, we have used finite-temperature 
ED/DMFT to investigate correlation effects in t2 g 
three-band transition metal oxides, such as Ca2Ru04, 50 
m x Go0 2 r^- LaTi0 3 , 52 and V 2 3 . 53 It was shown 
that, in these systems, accurate projections of the lattice 
Green's function onto a finite cluster consisting of three 
impurity levels and six bath levels (2 bath levels per 
t 2g orbital) can be achieved, yielding an overall cluster 
size n s — 9. Since the different baths are indirectly 
coupled via the 3d interorbital Coulomb and exchange 
interactions, the spacing between excitation energies 
is rather small, so that finite-size effects are greatly 
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FIG. 1: (Color online) Density of states components of 
LaFeAsO for tight-binding band structure derived in Ref. 24 . 
(a) t2 g components p X z, y z and p x 2_ y 2, (b) e g components p xy 
and p z 2 . 



bands due to Coulomb interactions. The local lattice 
Green's function is given by 

G mn (iuj n ) = V] (iu n + \i - H(k) - E(iu> n )) (2) 

— V / ran 

k 

where u) n = (2ri + l)nT are Matsubara frequencies and 
p, is the chemical potential. Since we consider paramag- 
netic systems, the spin index of G and E is omitted. As 
a result of the symmetry properties of H mn (k), the den- 
sity of states matrix is diagonal: p mn (uj) — 5 mn p m (ui). 
Local Coulomb interactions preserve this symmetry, so 
that G mn and E mn are also diagonal. We denote these 
components by G m (iu) n ) and E m (iw„). We point out 
that, because of the non-diagonal nature of H(k), each 
G m component is influenced by all components E m . 

For the purpose of the quantum impurity calculation 
within DMFT it is necessary to first remove the self- 
energy from the central site. This step yields the impu- 
rity Green's function 



Go,m(*w n ) — [G m (iuj n ) + E TO 



(3) 



Within ED/DMFT the lattice impurity Green's function 
Go is approximated in terms of an Anderson impurity 
model for a cluster consisting of impurity levels e m =i...5 
and bath levels £fc=6...i5: which are coupled via hopping 
matrix elements V m k- Thus, Go tTn (ioj n ) 
where 



diminished. Here, we generalize this approach to five 
impurity orbitals, each coupled to two bath levels, i.e., 
n s = 15. 

Denoting the tight-binding matrix elements by 
H mn (\s.) the interacting Hamiltonian is given by: 



H 



-^mri(k)cv| lkcr C„ker 

mnker 

"F ^ U n imf n imj. "i" ^ (U — JS aa i)ni m(7 ni m i n 
irn irn<m' era' 

— ^ J [cf m ^Ci m lcl m , ^Cim't + c tm J [ C tmi Cim ' / t Cim ' i\ 
ira^m' 

(1) 

where c^ CT are annihilation (creation) operators for elec- 
trons on site i in orbital m with spin a and riima = 

c tma c imcr- <^ma are tne corresponding Fourier compo- 
nents. The intra-orbital and inter-orbital Coulomb ener- 
gies are denoted by U and U'. The exchange integral is 
J, where U' = U — 2 J because of rotational invariance. 
Spin-flip and pair-exchange terms are denoted explicitly 
by J'. In the case of isotropic Hund exchange, one has 
J' = J. In the case of Ising-like exchange these terms are 
neglected, i.e., J' = 0. 

The aim of the five-orbital single-site DMFT calcu- 
lation is to derive the local self-energy matrix E mn (u;) 
which describes the modification of the single-particle 
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Since G§ m is diagonal in orbital indices, each impu- 
rity level couples to its own bath containing two lev- 
els: Orbital 1 couples to bath levels 6, 7, orbital 2 to 
bath levels 8,9, etc. Each of the four independent func- 
tions GQ m (iuj n ) therefore involves 5 adjustable parame- 
ters: one impurity level, two bath levels, and two hop- 
ping elements. These parameters are found by using a 
standard minimization procedure. The quality of these 
fits using 5 parameters is very good, as shown in several 
previous works i 53 ' 54 

As a result of the ED quantum impurity calculation 
one obtains the finite temperature cluster Green's func- 
tion which is also diagonal: G^(iw„). In analogy to 
Eq. §5§ the cluster self-energy is given by 



G 0,m( iuJ n)' 



(5) 



The important physical assumption within DMFT is now 
that this cluster self-energy provides an adequate repre- 
sentation of the self-energy of the extended solid, i.e., 
E^(io;„) s» E TO (io; n ), which is then used in Eq. (0) to 
derive the lattice Green's function in the next iteration 
step. Further details can be found in Ref.—. 

Since the cluster Hamiltonian is extremely sparse (typ- 
ically only 20 to 30 off-diagonal elements per row), the 
quantum impurity calculation is conveniently carried out 
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FIG. 2: (Color online) Fe 3d orbital occupancies (per spin) as 
functions of Coulomb interaction for n = 6, with J = [7/4 for 
U < 3 eV and J = 0.75 eV for U > 3 eV; Hund exchange. 



by using the Arnoldi algorithm^ The largest spin sector 
corresponds to (n^,ni) = (7,8), giving matrix dimension 
(15!/(7!8!)) 2 = 6435 2 . To reduce storage requirements, 
we have rewritten our multi-orbital ED/DMFT code so 
that large basis vectors of size 2 2 ™ s are avoided by keeping 
only vectors of size 2™ s . Moreover, the Arnoldi scheme is 
readily parallelized. Thus, using 32 processors the largest 
Hamiltonian subblock requires less than 1 GB storage. 
Since the spacing between excited states is very small, 
at finite temperatures a large number of states may con- 
tribute to the cluster Green's function. To reduce compu- 
tational time in this first five-band ED study, we perform 
the DMFT calculation at T = 0.01 eV, but retain only 
the lowest few states, making sure that ground state de- 
generacies are properly treated. Using 32 processors, one 
iteration then takes of the order of one to four hours. 



III. RESULTS AND DISCUSSION 
A. Undoped LaFeAsO 

An important consequence of local Coulomb interac- 
tions is the rearrangement of electrons among subbands. 
Fig. [2] shows the variation of the Fe 3d orbital occupan- 
cies with U at n = 6 total occupancy The precise values 
of U and J for LaFeAsO depend sensitively on the basis 
functions used for the 3d Hamiltonian . 18 ' 20 ' 21 ' 32 In the 
range U < 3 eV we chose J = U/4. For illustrative pur- 
poses we also show results for larger U. In this range J 
is kept constant at 0.75 eV, in order to avoid unrealisti- 
cally large Hund parameters. A charge flow from d z 2 to 
d x 2_ y 2 is seen to take place, thereby reducing the orbital 
polarization of the uncorrelated bands. The occupancies 



of the dL 



orbitals are less strongly affected by cor- 



relations. The results shown are for full Hund coupling. 
Ising exchange yields a similar charge rearrangement pre- 
dominantly between d z 2 and d x 2_ y 2 orbitals, with only 
slightly larger modifications of the d xz>yz>xy occupancies 
than seen in Fig. [5] Near U = 2.5 eV, all 3d occupan- 
cies are within about 10 % of the average occupancy 0.6. 
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FIG. 3: (Color online) Imagimary part of self-energy compo- 
nents E m (iu n ) as function of Matsubara frequency for several 
Coulomb energies, (a) and (b): Hund coupling; (c) Ising ex- 
change. The color coding is defined in panel (b). 



The origin of the unusual reduced d z 2 <x 2_ y 2 orbital po- 
larization is the complex bonding- antibonding shape of 
the density of states which yields n z 2 > n x 2_ y 2 although 
E z 2 > E x 2_ y 2. A correlation induced reduction of orbital 
polarization is also found in Na K Co02 which exhibits a 
similar pseudogap in the density of states as a result of 
the strong 3d— 2p hybridization in the planar geometry.— 
To illustrate the effect of Coulomb correlations on 
the Fe 3d bands in more detail, we show in Fig. [3] the 
self-energy components as functions of Matsubara fre- 
quency. In the case of Hund coupling, the initial slope 
of ImE m (iw„) ~ ijj n increases up to about U = 3 eV, 
until the quasi-particle weights are reduced to about 
0.15 . . . 0.3. Beyond this Coulomb energy, the self-energy 
components exhibit a finite onset of -0.3... -0.6 eV, in- 
dicating that states at the Fermi level acquire a finite 
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FIG. 4: (Color online) Quasiparticle weights Z m of Fe 3d 
orbitals as functions of Coulomb interaction for n = 6; Hund 
exchange. 
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FIG. 5: (Color online) Orbital dependent low-frequency scat- 
tering rates -y m as functions of Coulomb interaction for Hund 
and Ising exchange. 



lifetime. Since this onset is much larger than what is 
expected due to finite temperature, it implies a break- 
down of Fermi-liquid behavior. The loss of coherence is 
strongest for d xy and weakest for d z 2 . Ising exchange also 
gives rise non-Fermi-liquid behavior, except that the on- 
set occurs at about U = 2 eV, i.e., at considerably lower 
Coulomb energy than for Hund exchange. 

Fig. @] shows the orbital dependent quasi-particle 
weights Z m = 1/[1 — ImcE m (iuj) / dui^^o] as functions of 
Coulomb energy for Hund exchange. Despite the differ- 
ent 3d orbital occupancies and different ti g and e g band 
widths, all five bands are seen to exhibit a similar reduc- 
tion of Z m with increasing U. Ising exchange yields a 
slightly steeper decrease of Z m up to about U — 2 eV, 
beyond which all self-energy components show a finite 
onset. 

There exists strong experimental evidence that Fe 
pnictides exhibit an orbital dependent effective mass en- 
hancement, and an concomitant narrowing of quasiparti- 
cle bands, of about a factor of 2 to 3i 31 ' 34 i 56 ' 57 According 
to the results shown in Fig. [4] these experimental find- 
ings are compatible with the present five-orbital picture 
if U w 2.0 .. . 2.5 eV, with J w U/4 and full Hund's 
rule coupling. The system would then be just below the 
Fermi-liquid to non-Fermi-liquid phase boundary. For 
Ising exchange, on the other hand, these Coulomb and 
exchange parameters would imply non-Fermi-liquid be- 
havior. The above small values of U evidently reflect the 
fact that, in a purely 3d electron description, screening 
via As and O p states greatly reduces the Fe 3d Coulomb 
interaction. Correspondingly larger U are appropriate 
within a more refined dp description i 18 i 21 i 30 i 32 i 33 

Figure [5] shows the variation of the low- frequency scat- 
tering rates 7 m = — ImS m (io; — > 0) with Coulomb en- 
ergy. For Hund as well as Ising exchange, the onset of 
non-Fermi-liquid occurs for all bands at the same critical 
U. Whereas Ising coupling leads to a sudden rise near 
U 2 eV, Hund exchange shows a smoother increase 
near Pk3 eV. These results are consistent with previous 
ones for a degenerate two-band model^ which showed 
that Ising coupling yields a lower critical Coulomb energy 



and a change from a continuous to a first-order transition. 
The damping rates for Ising coupling near U = 2.5 eV are 
about 0.4 . . . 0.7 eV, in qualitative agreement with quan- 
tum Monte Carlo (QMC) results in Refs^S 2 .. For Hund 
coupling near U — 3.5 eV they are of similar magnitude. 
They continue to rise at larger U and indicate increasing 
orbital differentiation, with larger damping for d xy , X z,yz 
than for d x 2_ y 2 z 2. 

The transition from coherent to incoherent metallic 
behavior should also manifest itself in the temperature 
variation of the Fe 3c? self-energy. This has been studied 
recently for LaFeAsOi-^Fa; at x = 0.1 within continuous- 
time QMC DMFT for full Hund exchanged The co- 
herence temperature was shown to be a highly sensitive 
function of J, becoming extremely small for J w 0.7 eV, 
which is close to the value assumed here (J = 0.75 eV 
for U > 3 eV). 

To explore the origin of the Fermi-liquid to non-Fermi- 
liquid transition we have evaluated the spin-spin corre- 
lation function C mz (r) = (5 mz (r)5 mz (0)), where r de- 
notes the imaginary time. Fig. [6] shows these orbital de- 
pendent functions for several Coulomb energies and Hund 
coupling. At low values of U, the C mz (r) decay to zero, 
as expected for a Fermi liquid. The orbital components 
of the spin susceptibility 

/ dr (S mz {T)S mz (0)) (6) 
Jo 

are then independent of temperature, indicating Pauli 
behavior. With increasing U, the decay becomes less 
rapid and finite values are approached at large r (for r <C 
/3 = 1/T), demonstrating the formation of local moments 
S m simultaneously in all subbands. The susceptibility 
components are then proportional to j3, so that x-m ~ 
S m {S m + as expected for Curie- Weiss behavior. 

As shown in Fig. for Ising coupling the formation of 
frozen moments sets in at much lower Coulomb energies. 

A similar spin-freezing transition was recently found by 
Werner et al^- for a fully degenerate three-band model 
near n = 2 occupancy. Using continuous-time QMC 
as impurity solver, the paramagnetic U — n phase di- 
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FIG. 6: (Color online) Variation of orbital dependent spin- 
spin correlation functions with imaginary time r, for three 
Coulomb energies and Hund coupling, (a) U — 1.75 eV, (b) 
U = 2.5 eV, (c) U = 3.25 eV. 



agram was shown to exhibit Fermi-liquid properties at 
small U. For increasing U and n > 1.5, an incoher- 
ent metallic phase with local moments appears, which is 
then replaced by a Mott insulating phases at integer oc- 
cupancies n = 2 and n = 3. Beyond the critical value of 
U, the low-frequency limit of the self-energy exhibits a 
finite onset of similar magnitude as shown here in Fig. [5] 
In the present five-band study, we find in addition that 
this transition changes approximately from continuous to 
first-order when Hund exchange is replaced by Ising-like 
coupling. 

To illustrate the correlation-induced transfer of spec- 
tral weight in LaFeAsO we show in Figs. EJa) and (b) 
the 3d spectral distributions for two Coulomb energies 
and Hund coupling. For simplicity we plot here the ED 
cluster spectra since they do not require analytic continu- 
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FIG. 7: (Color online) Variation of orbital dependent spin- 
spin correlation functions with imaginary time r, for two 
Coulomb energies and Ising exchange, (a) U = 1.5 eV, (b) 
U = 3.0 eV; J = U/4. 



ation to real frequencies. The main effect at U = 1.75 eV 
is the band narrowing both below and above Ep. In ad- 
dition, spectral weight is shifted below the bottom of the 
3d bands, indicating the formation of weak lower Hub- 
bard bands. At U — 3 eV, spectral weight in the oc- 
cupied part of the 3d bands is greatly reduced and the 
Hubbard bands are much more prominent. Since there 
is little experimental evidence for any significant lower 
Hubbard bands, the above results imply that U within 
the present 5-band description should be approximately 
2.0 . . . 2.5 eV, where the precise value depends on the 
magnitude of J. The unoccupied states are also shifted 
closer to Ep , but there is no evidence of any upper Hub- 
bard peaks. Note that the lower Hubbard bands are very 
broad. This is related to large uncorrelated band width 
and to the multiplet structure induced by Hund's rule 
coupling. The spectra for Ising exchange shown in panel 
(c) are qualitatively similar. Many small differences arise 
because of the absence of spin flip and pair hopping ex- 
change processes. The different multiplet structures as- 
sociated with Hund and Ising coupling will be considered 
in more detail elsewhere. 

The above results demonstrate the importance of a 
proper treatment of exchange interactions. For instance, 
if because of pd screening realistic values of U and J for 
LaFeAsO are approximately 2.0 . . . 2.5 eV and J « U/4, 
respectively, full Hund coupling suggests that the system 
is moderately strongly correlated, with 3d effective mass 
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FIG. 8: (Color online) Spectral distributions of Fe 3d sub- 
bands for (a) U — 1.75 eV, (b) U = 3 eV; Hund coupling; 
(c) U = 2.5 eV: Ising exchange; with J = U/4 and 0.05 eV 
broadening. 



enhancements of the order of m* /m = 1/Z m « 2 ... 3. In 
contrast, if spin- flip and pair-exchange processes are ig- 
nored (as is usually done in Hirsch-Fye QMC calculations 
to avoid sign problems, e.g., in RefSi 20 ! 21 ' 33 ), the same 
Coulomb and exchange parameters suggest that system 
has crossed the boundary towards non-Fermi-liquid be- 
havior, with strongly reduced lifetimes of electronic states 
close to Ep. It would be interesting to inquire whether 
a dp formulation of i?(k) yields a similar qualitative dif- 
ference between Hund and Ising exchange. 

At sufficiently large Coulomb energies, n = 6 integer 
occupancy should eventually lead to a Mott insulating 
phase. We have increased U up to 6 eV while keeping J = 
0.75 eV constant. Both for Hund and Ising coupling, the 
system evolves towards an orbital selective phase, where 
the d xz .y Z . X y subset is either in or close to an insulating 



phase and the d x 2_ y 2 z 2 subset remains in the strongly 
incoherent metallic stated Thus, for realistic Coulomb 
and exchange we conclude that in the present five-orbital 
ED /DMFT description the system is far below the n = 6 
Mott insulating region. Orbital selective phases in five- 
band systems at n = 6 occupancy, with subbands split 
via a crystal field, were also found in Refs i 27 ' 59 

We close this subsection by pointing out that several 
papers have recently discussed the role of coexisting itin- 
erant and localized electrons in iron pnictides^ ' 10 i 26 ~ 28 i 36 
A simple system which exhibits this kind of coexis- 
tence is the half-filled two-band model where the nar- 
row band is Mott localized while the wide band is still 
metalliCi 46 i 60 ~ — Coulomb interactions between these two 
types of electrons give rise to bad-metallic behavior in 
the itinerant band, with a finite scattering rate and a 
pseudogap at Ep in the case of Ising exchange,— and 
marginal Fermi-liquid behavior in the case of full Hund 
coupling] 61 ' 62 As we have seen above, the important fea- 
ture of FeAsLaO at n = 6 occupancy is that the spin 
freezing transition occurs simultaneously in all five d 
bands, far below the Mott transition. Thus, in the range 
of Coulomb and exchange energies of interest, all sub- 
bands have similar occupancies, effective masses, and 
low-frequency scattering rates. The subbands therefore 
do not split into itinerant and localized subsets. Only at 
much larger U, the d XZtVZ>xy and d x 2_ y 2 z 2 bands gradu- 
ally tend towards 0.5 and 0.75 occupancy. 



B. Doped LaFeAsO: correspondence between 
pnictides and cuprates 

The results discussed above are for paramagnetic, un- 
doped LaFeAsO, with Fe 3d occupancy n — 6. To il- 
lustrate the effect of electron and hole doping we show 
in Fig. IH1 the orbital occupancies as functions of chem- 
ical potential at fixed U — 2.5 eV. In order to eluci- 
date the doping variation of many-body properties, the 
one-electron Hamiltonian is kept unchanged. The doping 
range extends from one hole to one electron. (At half- 
filling n = 5, the chemical potential is fj, = 4.5U — 10 J = 
5 eV.) Evidently, the degree of orbital polarization de- 
pends strongly on the total occupancy. For n = 5, all 
d bands become half-filled and some orbital components 
ImE m (iw ra ) are porportional to l/uj n . The spectral dis- 
tributions reveal that the system then is a Mott insulator 
where all d bands are split into lower and upper Hubbard 
peaks. 

In the case of electron doping, all orbital occupancies 
increase in a similar fashion. At n = 7, the self-energies 
(not shown) indicate that the system is a weakly corre- 
lated Fermi liquid. Even for Ising exchange, the quasi- 
particle weights are in the range Z m = 0.26... 0.56. 
These properties differ strikingly from those discussed 
above at n = 6, where under the same interaction condi- 
tions the system is much closer to bad-metallic behavior. 
Thus, the Fermi-liquid to non-Fermi-liquid phase bound- 
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FIG. 9: (Color online) Fe 3d orbital occupancies as functions 
of chemical potential for U = 2.5 eV with J = 0.625 eV and 
Hund coupling. The arrows indicate the values of fi associated 
with integer total occupancies. 



ary for n — 7 is shifted to larger values of U. These results 
suggest a fundamental asymmetry of LaFeAsO with re- 
spect to doping. For hole doping, bad-metallic behavior 
should increase, while electron doping reinforces Fermi- 
liquid properties. This behavior is consistent with re- 
sults for the three-band model^ where the Fermi-liquid 
to non-Fermi-liquid phase boundary was also found to 
shift to larger U as the occupancy moves farther away 
from half-filling. 

On the basis of the above results we obtain the para- 
magnetic U — n phase diagram shown in Fig. 1101 where, 
for U = 2.5 eV and J = f//4 Hund exchange, undoped 
LaFeAsO lies just below the Fermi-liquid / non-Fermi- 
liquid phase boundary. Smaller U and larger J would 
move this point farther below this phase boundary. In 
the limit of one-hole doping, the system is a Mott insula- 
tor, whereas, for electron doping, Fermi-liquid properties 
dominate. For Ising exchange, the phase boundary is 
shifted to roughly 1 eV lower U values, so that U would 



4 



> 









! / 1 

Hunri/ y 




Mott 


NFL 


■r Ising ■ ■ 








LaFeAsO i 






FL 


i . i 



5 6 7 

n 

FIG. 10: (Color online) Schematic phase diagram for doped 
LaFeAsO. Red curves: boundaries between Fermi-liquid (FL) 
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change. At half-filling (n = 5) a Mott insulating phase exists 
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7 m for LaFeAsO as functions of doping relative to half-filling, 
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nents of scattering rate for single-band two-dimensional Hub- 
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nearest neighbor hopping t = 0.25 eV and t' = —0.075 eV, 
respectively, U = 2.5 eV, T = 0.01 eV), 5 - The vertical bars 
denote the approximate doping concentrations of the Fermi- 
liquid to non-Fermi-liquid transition. 



have to be less than ~ 2 eV to preserve Fermi-liquid be- 
havior. 

Of course, in real LaFeAsO, this phase diagram should 
be more complicated because of the inevitable modi- 
fication of the one-electron properties with dopin g 15 i 63 
and because of the antiferromagnetic phased observed at 
n = 6. Nevertheless, the paramagnetic limit permits to 
draw an interesting analogy between the present multi- 
band iron pnictide and the single-band two-dimensional 
Hubbard model that is frequently used to investigate 
Coulomb correlations in high-T c cuprates. 

Fig- [TIT a) shows the orbital components of the 
low-frequency scattering rate 7 m « — ImS m (iwo) for 
LaFeAsO as a function of electron doping relative to half- 
filling, i.e., S = n/5 — 1, where n = 6 for the undoped 
material. As discussed above, in the limit of one-hole 
doping (n — 5), the system is a Mott insulator, while 
for S > 0.2 it becomes a weakly correlated Fermi liquid. 
In the intermediate region, for S < 0.2, the scattering 
rate increases sharply so that the system is dominated 
by non-Fermi-liquid properties due to the formation of 
frozen moments. Thus, close to n = 6 the system is near 
the Fermi-liquid to non-Fermi-liquid phase boundary. 

As shown in panel (b) , this behavior is remarkably sim- 
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ilar to the one found within cluster ED/DMFT for hole 
doping in the two-dimensional Hubbard models Anal- 
ogous results have been obtained by several groups.— ~— 
This model yields a Mott insulator at half-filling and ex- 
hibits a non-Fermi-liquid pseudogap phase up to a criti- 
cal doping S c w 0.15 . . . 0.2. Ordinary metallic behavior 
is restored in the overdoped region, S > 6 C . For electron 
doping the results are similar, except for a smaller critical 
doping which marks the onset of bad-metallic behavior. 
In view of the analogy seen in Fig. [TTJ the Mott phase 
that is relevant for LaFeAsO is not the one that should 
eventually appear at n = 6 for large U, but the one that 
exists at realistic values of U and J at n = 5 occupancy. 

The results discussed above demonstrate that the Mott 
transition in multi-orbital and multi-site Hubbard mod- 
els is far more complex than the paramagnetic metal to 
insulator transition obtained in single-band, single-site 
DMFT calculations. The presence of interactions chan- 
nels involving orbitals or sites not only affects the over- 
all magnitude of the critical Coulomb energy, but gives 
rise to a much richer phase diagram. In particular, as 
we have seen here for LaFeASO, a frozen-moment, non- 
Fermi-liquid phase emerges between the Fermi-liquid and 
insulating regions. Conceptually, this non-Fermi-liquid 
behavior is closely related to the pseudogap phase in hole- 
doped cuprates which arises from short-range Coulomb 
correlations. 



IV. SUMMARY 

Multi-band ED/DMFT has been used to investigate 
the effect of correlations in the iron pnictide LaFeAsO. 
Starting from an accurate five-band tight-binding single- 
particle Hamiltonian, the many-body properties are eval- 
uated by extending single-site ED/DMFT to five orbitals, 
each hybridizing with two bath levels. This scheme 
is particularly useful at very low temperatures, large 
Coulomb energies and for fully rotationally invariant 



Hund's rule coupling. It is shown that correlation ef- 
fects in LaFeAsO give rise to a paramagnetic transition 
from a Fermi-liquid phase to a non-Fermi-liquid phase at 
a critical Coulomb energy (pd screened) of about 3 eV. 
This transition appears to be continuous and is caused by 
the formation of Fe 3d local moments. Below this transi- 
tion, the quasi-particle weight of all d orbitals is strongly 
reduced and orbital polarization is less pronounced than 
in the uncorrelated density functional band structure. A 
Mott insulating phase does not appear below [7 = 6 eV. 
For Ising exchange the Fermi-liquid to non-Fermi-liquid 
transition also exists, but the critical value of U is shifted 
to about 2 eV. Moreover, the transition appears to be 
first order rather than continuous. 

Despite the multi-band nature of LaFeAsO and the 
important role of exchange interactions, this system ex- 
hibits an interesting relationship to the single-band two- 
dimensional Hubbard model if the doping concentration 
is defined with respect to the half-filled 3d shell. Within 
this picture, iron pnictides and cuprates exhibit the same 
sequence of paramagnetic phases with increasing doping: 
from Mott insulator to bad metal to Fermi liquid. Thus, 
fluctuations between orbitals or lattice sites lead to sim- 
ilar physics. According to the available experimental ev- 
idence, LaFeAsO is located on the weakly to moderately 
correlated side of the Fermi-liquid to non-Fermi-liquid 
phase boundary. Thus, iron pnictide materials appear 
to have two parent compounds: the anti-ferromagnetic 
semi- metal at n — 6 and the Mott insulator at n = 5. 
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